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In this paper we construct an effective field theory for a condensate of cold Fermi atoms whose 
scattering is controlled by a narrow Feshbach resonance. We show how, from first principles, it 
permits a hydrodynamic description of the BEC-BCS crossover from which the equation of state, 
intimately related to the speed of sound, can be derived. Specifically, we stress the non-universal 
behaviour of the equation of state at the unitary limit of infinite scattering length that arises when 
either, or both, of the range of the inter-atomic force and the scale of the molecular field become 
large. 



I. INTRODUCTION: SCATTERING WITH 
FESHBACH RESONANCES 

Cold alkali atoms whose scattering is controlled by a 
Feshbach resonance can form diatomic molecules with 
tunable binding energy on applying an external mag- 
netic field Weak fermionic pairing gives a BCS theory 
of Cooper pairs, whereas strong fermionic pairing gives 
a BEC theory of diatomic molecules. The transition is 
characterised by a crossover in which, most simply, the 
s-wave scattering length as diverges as it changes sign 

a. 

A considerable theoretical and experimental effort has 
been expended on understanding such macroscopic quan- 
tum systems. We wish to show that, in a formalism in 
which the tuning of the system by the external field is 
explicit because of the narrowness of the resonance, it is 
straightforward to derive the semiclassical attributes of 
the condensates ( speed of sound, hydrodynamics, equa- 
tion of state) analytically. Our starting point is the exem- 
plary 'two-channel' microscopic action (in units in which 
h = 1) 



S 



dtd 3 - 



u 

+ V^(x) i/>l(x) i/>i(x) ^t(x) 
V 2 



V 2 

i d t + h n 

2m 



g [<j>*{x) ^{x) V T (aO + fa) ^(x) VK*)] (1) 



for cold fermi fields ip a with spin label a = (1,1), which 
possess a narrow bound-state (Feshbach) resonance with 
tunable binding energy v, represented by a diatomic field 
(f> with mass M = 2m [H, [j, 0, Q . In addition there is 
a simple s-wave attractive contact interaction U > 0, 
regularised at a momentum scale A. Standard renor- 
malisation methods 0] allow us to remove A from the 



formalism. In what follows we assume that such renor- 
malisation has been made. The effective atomic coupling 
strength (at zero external energy-momentum) is 



U eff = U + g 2 /{v - 2/i) = -k F a s /N 



(2) 



comprising the contact term plus the effect of the res- 
onance, where Nq is the density of states at the Fermi 
surface. Tuning v — 2fi to zero by the application of an 
external magnetic field sends \as\ — > oo. This is the so- 
called 'unitary' regime. 

This work follows on from that of an earlier paper 0] , 
henceforth referred to as I, from which we recreate the 
condensate effective action © in the next section. 

We shall show, in greater detail than in I, that the 
(T = 0) condensate described by |T]) can be understood, 
in the hydrodynamical approximation, as two coupled flu- 
ids, built from the fcrmion-paired atomic and molecular 
subsystems. This is sufficient to determine the equation 
of state (EOS) and speed of sound exactly in the mean- 
field approximation. There are simplifications in that a 
single fluid dominates in a) the deep BEC regime b) the 
deep BCS regime and c) the unitary regime. When a 
single fluid description is appropriate the EOS can show 
several different allometric behaviors p oc p l+1 depend- 
ing on the magnitudes of the effective range of the inter- 
atomic force and the length scale of the molecular field. 

In particular, when either is large, the single fluid be- 
haviour in the unitary regime does not show the canonical 
value 7 = 2/3 and we do not have a conformal field theory 
when | ag | — + oo. There has been much activity recently 
on taking advantage of the dualities of Anti de-Sitter gen- 
eral relativity and conformal field theory (AdS/CFT) Q 
to use classical black hole physics to describe cold Fermi 
atoms in the unitary regime in those cases when the the- 
ory is conformally invariant Q , but our models show the 
limitations of this analysis. 

Possible experimental tests are considered. 
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II. EFFECTIVE ACTIONS 



rcscalc <J|A| to 5\A\ = ne, such that L e ff takes the form 



Introducing auxiliary fermion-paired bosonic fields 
A(x) = Uijji{x)^{x), A*(x) = Uipl(x)^(x) renders 5 
quadratic in the fermi fields. Integrating them out Q 
enables us to write S in the non-local form 



-iTrlng- 1 + j dt dx( - ^|A| 



4>*{x) 



2M 



* &t + 7777 + 2/i — V 

in which Q^ 1 is the inverse Nambu Green function 

q-1 _( idt-e A(x) 
1 A*(x) idt+e 



4>{x)\, (3) 



(4) 



where A(x) = A(x) — 3 represents the two- 

component combined condensate (and £ = — V 2 /2?n— /.t). 

In this paper we restrict ourselves to the mean-field 
approximation, the general solution to SSnl = 0, valid if 
is a sufficiently narrow resonance [H, . Our approach 
is a generalisation of that in [ll|, discussed by several 
authors from [l2j onwards (see Schakel (HI). 

If we write A(x) = \A(x)\ e i0 ^ and cf>(x) = 
— \4>(x)\ e 18 *^ the combined condensate amplitude and 
phase of A{x) = |A(x)| e l8 ^ x ^ are then determined. 
The action possesses a U(l) invariance under 0a — > 
0A + const., — ► fl^ + const., which is spontaneously 
broken: (5 S/vl = permits spacetime constant gap solu- 
tions |A(x)| = |A | ^ and \<j>(x)\ = \<j) \ ^ (whereby 
(|A(ir)| = | Ao| ^ 0) and a Goldstone boson, the (gaplcss) 
phonon. |Ao| determines the density of states at the 
Fermi surface as N = / d 3 p/(27r) 3 (| A | 2 /2£^), where 



El = et 



|A | 



s p = pV2m-/i. If |A | = |A | + |-^o|, 
then |A |/|A | = U/U eS , \c/> \/\A \ = g/[{v - 2p) U eS ]. 
In addition, the system possess a gapped ('Higgs') mode. 

To determine the EOS it is sufficient to consider just 
the fluctuations around the gap configurations and per- 
turb in the small fluctuations in the scalar condensate 
densities [O, E31 ^l A l = l A l ~ l A ol and 5\tj>\ = \<j>\ - |0 O | 
and their derivatives. We perform a Galilean invariant 
long wavelength, low-frequency expansion in space and 
time derivatives to give 



S e ff = I dtdxL 



J eff 



(5) 



in terms of the local effective Lagrangian density L e ff 
with elliptic equations of motion. Although a and 00 
are not small, we assume that <5|A|, 8\4>\, and (0a — &<p) 2 
are of the same order. 

Because Q^ 1 is defined in terms of A it is natural to 
express the local Lagrangian density L e f f in terms of the 
phase angles A , 00 and the Galilean scalar fluctuations 
<5|A| = |A| — |Aq| and 6\</>\. In fact, it is convenient to 



Leff = -\ P BG{6^)- l -tf{6 K -e^ 

-i(^ + 2« £ )G(0 A , e ) + ^G 2 (0 A ,6) 
+ \nX 2 {e,9 K )-\M 2 e 2 + ^5^ 

£ + (2,-^)^|, (6) 

in terms of the Galilean scalar combinations G(0) = 
+ (V0) 2 /4m, G(0,e) = + (V0) 2 /4m + (Ve) 2 /4m, 
X(e, 0) = e + V0.Ve/2m and A - 0$. We have chosen 
k so that the dimcnsionlcss e has the same coefficients as 
A in its spatial derivatives. On taking g identically zero 
in © we recover the Lagrangian of [ll| ■ 

The fermion number density arising from the gap 
equations is p° = p F + p B , where p F = 
J d 3 p/(2ir) 3 [1 — Ep/Ep] is the explicit fermion density, 
and p° B = 2 1 0o | 2 is due to molecules (two fermions per 
molecule). We shall introduce the other coefficients, 
which arc simple momenta integrals, as and when they 
are needed. It is the form of © rather than the detail 
that concerns us at the moment. 



III. EQUATIONS OF STATE 



The definition of a narrow resonance is that 



7r 



v / r7^ : « i, 



(7) 



where F is the resonance width and ep = hp /2m is the 
typical atomic kinetic energy. Unless stated otherwise, 
we assume narrow resonances, for which the mean field 
approximation can be justified [icj |. Since the Fermi mo- 
mentum kF increases as the density p increases we can, 
in principle, make even broad resonances narrow by in- 
creasing the density. In more detail [IcI ]. 



m 2 {2p-v) 2 U 2 ff 



(37r 8 )V3 



p l/3 ff 2 



(8) 



Note that j r oc g 2 when \as\ — > oo, irrespective of U. 
However, j r increases as we move into the deep BEC 
and BCS regimes, when \v\ is large, making the narrow 
resonance approximation less valid. With this qualifica- 
tion, systems with narrow resonances include e Li with 
7 r w 0.2. On the other hand, when 7 r ^> 1 and the 
narrow resonance approximation breaks down, the model 
effectively becomes a one-channel model in its basic prop- 
erties. There are then strong similarities with the single- 
channel model in which g = identically, discussed by 
many authors, but for which we cite [r"0 ] in particular. 

The hydrodynamics of the system is encoded in 
9&(x) and 00 (cc) and will have a natural realisation as 
two coupled fluids. To proceed, we ignore the density 
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and velocity fluctuations p e = ~(rj/A)X(e : 9^) and 
v e = Ve/2m due to the condensate fluctuations e, in 
comparison to p° and v = V#^/2m, the condensate 
velocity. The inclusion of p e and v e in the two-fluid 
model would give small fluctuating short-range sources 
and sinks in the fluids. All that we need for the EOS is 
the hydrodynamic approximation, which coarse-grains 
by replacing them with their (zero) averages. 

The angular Euler-Lagrange (EL) equations are 
then 



0_ 

Of 
d_ 

at 



PF + V{p F v)-2n 2 (6 A -6 4> ) = 0, 
PB+^(pBu) + 2n 2 (0 A -9 4> ) = 0, (9) 



where u = V8<p/2m and pf = P° F — N G(9^,e) + 2cee, 
p B = 2\(f>\ 2 w p% + 4\(f> \5\(j)\. Putting these together 
gives 



^{Pb + Pf) + V.(p B u 



p F w) 



0, 



(10) 



the continuity equation for two coupled fluids, as given 
in I. 

Eq. ljlOp can be written in a more transparent form. 
Whereas the explicit fermion density in molecules ps in 
([TO]) appears in conjunction with the velocity u of the 
molecular component of the fluid, the generalised fermion 
pair density pp is coupled to the combined condensate 
velocity v, rather than v = V9&/2m. However, from 
the definition of A it follows that 9^ = b9& + (1 - 
where b = |A |/|A | = U/U e ff. In consequence, v 
bv + (1 — b)u, whereby (|10[) can be written as 



d_ 

at 



(pb + Pf) + V.(p B u + p F v) = 0, 



(11) 



where p F = bp F and ps = Pb + (1 — b)p F . That is, the 
effective molecular density ps describes point-particle 
bosons together with a cloud of fcrmionic Cooper pairs, 
which deplete the effective fermion pair density p F . Be- 
cause the fluids are coupled, the condensate moves as a 
single entity with velocity v. 

The Bernoulli equations from which the EOS follows 
are derived from these EL equations on substituting for 
the densities. Again neglecting p e and v e , the EL equa- 
tion for 9^ can be written as the simple Bernoulli equa- 
tion 



Shp + —m v 2 



= 0, 



(12) 



where Sh F = (p F — p F — 2ae)/2N is the specific enthalpy. 
After substituting e from its EL equation 

M 2 e + 2a G{9 K ) - (Ag/U)6\(f>\ «0 (13) 

in Sh Fl the enthalpy can be expressed in terms of the 
density fluctuations Sp F = p F — p F and 8ps = 4\<fio\6<f> 



as 

where 
K FF 



dp 

Sh F = — = K FF 5p F + K FB Sp B (14) 
Pf 



M 2 



K 



2(N M 2 + 4a 2 ) ' " FB U \ (j> Q \{N M 2 + 4a 2 )' 

(15) 

Complemcntarily, the EL equation for 9a, has the form 



-get 



raw. 



Sh f 



—mu — 



2"" 16mp0j 



V 2 5 Pi 



= 0, 



where hs permits the decomposition 



8p 



Sh,B = -fr = KbfSpf + KbbSpb- 
Pb 



As required, Kbf = Kfb and 



K B b = 



{v - 2/i) U e ff 
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No 



U 2U 2 \(t> \ 2 4a 2 + N M 2 ' 



(16) 



(17) 



(18) 



The \7 2 Sps oc V 2 5|</>| term in (TIB)) is just as we would 
expect from a theory of a pure bosonic gas. In the hydro- 
dynamic approximation such derivatives of d\(f>\ are also 
ignored comparatively, and the resulting equation 



rau + V 



5hi 



= 0, 



(19) 



is taken in conjunction with (|12[) in determining the EOS. 
It follows from (H]) and JT7} that 

6p F (p°pKpp- p%K BF ) = 5p B (p%K B B- P F Kp B ). (20) 

We have learned from (fT0|) and (jTTJ) that, however it may 
be partitioned, the total fermion density is p F +PB, whose 
fluctuation is Sp = Sp F + 5ps- Eqs. (fbl|) and (JTTJ) then 
collapse to give the EOS of the condensate as 



(21) 



dp _ p F p%{K FF K B B - K FB K BF ) 

dp (p F K FF + p° B K B B - P°b k bf - P f K F b 



where we remember that K F b = Kbf < 0. Eq. (|2ip is 
the key equation of this paper. 

As a very good check on our calculations, the time 
derivatives of the EL equations for 9a, and 9a, determine 
the dispersion relations of the modes. It follows directly 
that the speed of sound v is given as 



K FF K 



FF^BB 



K 2 
^BF 



P" 



m V K FF + K B b - 2K B f 



(22) 



Neither dp/dp nor v 2 depend on the coefficients rj, fl 2 
in ^ or the scaling parameter k and we use the re- 
sults of 0] to identify the coefficients in the Ks as 
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a = fd 3 p/(2TT) 3 (\A \e p /2E 3 p ) andM 2 
where 



2(2/U-P) > 0, 



d 3 p 

(2tt) 3 



El (p2/2m) 



< 0. 



(23) 



We assume that all the parameters in dp /dp and v 2 have 
been renormalised as in I. On substituting this and the 
Ks above into v 2 of (|22|) we do, indeed, recover the results 
of I. Unfortunately, although dp/ dp is given in terms of 
straightforward momenta integrals, in general it bears no 
simple relationship to v 2 . This is not surprising since the 
system comprises two coupled fluids. 

We conclude this section with a further comment 
on the role of the order parameter fluctuations e and 
6\4>\ (deviations from homogeneous gap solutions). The 
coarse-graining described above is not only sufficient to 
obtain the EOS, but obligatory within the approximation 
(J3J), that corresponds to taking only the lowest relevant 
field derivatives, in that the fluctuations are too short- 
range or too fast [llj to be properly accounted for by 
it. 

To go beyond the EOS requires an understanding of 
the fluctuations of the gapped 'Higgs' mode mentioned 
earlier, built from order parameter fluctuations. (This 
is trivially so in the limi ting case of g = for which 
the Higgs field is just e JTlj] . ) Equally, the dynamics 
of the Higgs mode cannot be obtained reliably from the 
approximation |(5J) and we need to incorporate higher- 
order field derivatives non-perturbatively to describe it 
(e.g. see flU}). 

For our purposes we can ignore the Higgs mode and 
we revert to our results of (|2ip and 



IV. SINGLE FLUID REGIMES 

With the qualifications above, there are circumstances 
in which we have a description in terms of a single fluid. 
Most simply, in the deep BCS regime (where ps ~ 0) and 
the deep BEC regime (where pp m 0), the system of (|10p 
behaves as a single fluid, of Cooper pairs or molecules 
respectively. A characteristic of these extremes is that 
the terms that are negligible in dp/ dp are also negligible 
in v 2 , whereby we get the simple result dp/ dp oc mv , as 
we shall now show. 



A. The BCS regime 



In the deep BCS regime p B = 2|0o| 2 vanishes, as does 
a because of particle- hole symmetry. As a result, Kbb 
becomes very large and 



dp 
dp 



ppKpp » -B- 
HF 2A 



m v 2 . 



(24) 



The EOS follows directly. At the Fermi surface (p 
the density of states Ao ex mp 1 / 3 whence 

dp 2/3 
— oc p oc p 1 
dp 



(25) 



or, equivalcntly, p oc p 1+1 where 7 = 2/3. There is no 
pressure from molecules in this limit. This is in agree- 
ment with the work of other authors (e.g. (l7l. [l8j) 



B. The BEC regime 

The deep BEC regime (large negative v) is charac- 
terised by small No, large a, and M 2 « 4/U. Un- 
like the BCS regime, Kpp dominates the Ks, so that 



PbKbb /m. Inspection shows that 



dp 
dp 



—m v 
2 



(26) 



since p F Kpp « p B \Ksp\- We have shown in I that in 
the BEC regime v 2 oc p, giving us the EOS 



p oc p 



1+7 



with 7 = 1. 



(27) 



This is the expected value of 7, corresponding to a molec- 
ular condensate with repulsive interaction [191 ] . 



C. The Unitary limit 

The extremes above permit relatively trivial single- 
fluid descriptions. More important is the 'unitarity' limit 
describing the central region v = 2p where \ag\ — > 00. 
This is described by a single fluid, not because ps or 
Pp vanishes, but because b = U/U e j f — » 0. This is es- 
sentially a one-channel system [201, represented by the 
single fluid pF = 0, pB = p in terms of the ps of 
We stress that the divergence of \as\ is not a signal of 
singular behaviour. The K's that define the behaviour 
of the system vary continuously as we pass from the BCS 
to BEC regimes through the unitarity regime, and vice- 
versa (e.g. see (flB)) ). for which pcrturbative methods 
such as ours are valid for narrow enough resonances (e.g. 
see Fig. 4 of [3). 

In addition to the scattering length as, the general 
model possesses two other important length scales, the 
effective range of the force ro where -f r ~ h/kp\ro\ as 
well as the length scale £ = (|4m (2p - v) Ueff/U])' 1 / 2 
for the molecular field. In general both ro = 0(g~ 2 ) 
and £ = 0((g 2 /C/) -1 / 2 ) are small near the unitary limit 
for broad resonances with large resonance coupling g, 
and the only effective dimensionless parameter is fcp|as|. 
On it diverging, the theory shows universal behaviour, 
with conformal invariance. In principle, AdS/CFT dual- 
ity then enables us to convert the difficult strong cou- 
pling calculations needed to describe the system into 
more tractable weak coupling boundary calculations for 
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classical black holes |8j . However, for narrow resonances 
| ro | becomes large and this universal behaviour is bro- 
ken. Further, for small g 2 jU we see that £ becomes large, 
again leading to the breakdown of universal behaviour. 

Nonetheless, in many regards there is still de facto uni- 
versality for realistic systems [see Figs. 30 of Chen et al. 
[HI) for many observables. To see deviations from non- 
universality the best way is to look for the deviation of 
the EOS from its canonical behaviour. 

After UV renormalisation, the gap equation at this 
limit becomes 







1 



u, 



eff 



d 3 p 



1 



1 



This fixes 



(2n) 3 [2E p 2(p 2 /2m) 



|A | = cfx, 



(28) 



(29) 



where c f=a 1.16. In turn, this relates the density p to p 
as 

2c 2 "> c l , n s3/2 



P = PB + PF- g2 r , ^ 



-^P 2 + ^rnp) 3 ' 2 



(30) 



where the first and second terms are ps and pp respec- 
tively and ci w 1.47. 

Although we have constructed the model with narrow 
resonances in mind the formalism permits extension to 
broad resonances [ll[ and we consider all possibilities. 
In the first instance, for broad resonances (7,. 3> 1), for 
which fc^|r | <C 1, p ~ pp pb, whereas for narrow 
resonances (7,, <C 1), for which ^fI^oI > 1, p w pB > 
pp. 

In determining the EOS wc stress that, in the unitary 
limit, the resonance width is independent of U. 



1. Broad resonances 

The canonical conformal symmetry in the unitarity 
regime arises when both fci?|ro| <C 1 and kp^ < 1 or, 
equivalently, we have large resonance coupling g, with 
g 2 > Up , p 2 /(2mp) 3 / 2 . For such values pfK ff <C 
PbKbbj p\Kbf\i and we recover the 'universal' behav- 
ior 



dp 2/3 
— oc n oc p 1 , 
dp 



(31) 



invariance is broken by kpt; 3> 1. In those circumstances 
PfK ff > pbKbb, p\K B f\, whereby 



d P 
— oc p oc p . 

dp 



(32) 



giving the EOS p oc p 1+1 with 7 = 0. 



2. Narrow resonances 

When the resonance is narrow the mean field approx- 
imation is robust [lo| . There is no contradiction in 
having strong self-interaction U (U >> pj (2mp) 3 l 2 ^> 
g 2 / p), when the conformal symmetry, already broken by 
^f|^o| 3> 1, is further broken with 



-1/2 



( P 



1/3 /g 2 \ -1/3 



> 1. (33) 



\ (2mp) J \ p 
In this case we find ppKpF 3> p\Kbf\ PbKbb, giving 



d P 
— oc p oc p . 
dp 



(34) 



That is, p oc p 1+7 where 7 = 0. This is the behaviour of 
a free Bose gas, expected in this limit (2l| . 

However, for relatively weak self-interaction 
(U,g 2 /p -C p/(2mp) 3 / 2 ), with smaller k F ^ than 
previously, we find pbKbb -C PfKff, p\K bf\, leading 
to 



dp 1/2 1/4 
— oc p ' oc p 1 . 

dp 

That is, p oc p 1+y where 7 = 1/4. 



(35) 



As we anticipated earlier, in all of the cases listed 
above it happens that dp/ dp oc mu 2 , and we can equally 
well read off the behaviour from that of v 2 . Away from 
these extremes we see, from (f3"0")) . that an allometric 
representation is not justified, although the effective 
exponent "f e ff = d{\ndp/dp)/d\np will interpolate 
between these values. 



or p oc p 1+1 with 7 = 2/3, as in the BCS regime. 

Although our approximation is not wholly reliable in 
that, for broad resonances, the renormalization of the 
molecular boson is expected to contribute sizable correc- 
tions to the EOS in such a strongly coupled regime [lj|, 
this is the correct value ljj , observed experimentally (sec 
Fig. 15 of [13). 

For broad resonances, p ~ kp/2m, whereby kp^ ~ 
(JJ pi g 2 ) 1 ! 2 . For relatively smaller g (Up ^> g 2 » 

p 2 /(2mp) 3 / 2 ), even though |ro | remains small, conformal Consider an elongated axially symmetric harmonic 



V. EXPERIMENTAL TESTS 

The best experimental test of our non-canonical 
equations of state is to determine the EOS exponent 7 
directly by observing the expansion of the condensate in 
elongated traps on removing the potential [TH, [l7| ■ This 
is governed by the hydrodynamical equations and to see 
how this carries over to our model we need to examine 
the single fluid nature of the unitary regime in more 
detail. 
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FIG. 1: The aspect ratio as a function of the time for the 
expansion of the condensates with different allometric be- 
haviours for its EOS when A = u)^_/u z = 0.1. The 'universal' 
behaviour 7 = 2/3 corresponds to strongly coupled broad 
resonances. For narrow resonances in the unitary limit 7 in- 
terpolates between 7 = and 7 = 1/4 and 7 = 1 arises in the 
BEC regime. 



trap with frequencies w x = oj y = uj±_ and lj z (u>± > oj z ). 
When there is simple allometric behaviour 5h oc p 7 , the 
density can take the scaling form 



(36) 



where b x = b y = b±. The continuity and Bernoulli equa- 
tions then reduce to [13 



l >> II, b l 



(37) 



The aspect ratio, the relative transverse to longitudinal 
expansion of the condensate, 



R±(t) = R±{0) b ± (t) 
R z (t) R z (0) b z (t) 



(38) 



can be measured where R± and R z are radial and axial 
radii respectively. In Fig. (Q} we show this ratio for the 
values of 7 derived above (in the unitarity and extreme 
BEC and BCS regimes) for A = ljj_/ui z = 0.1. It is 
evident that the hydrodynamical effect is greater in 
the direction of larger density gradients, as anticipated. 
For a larger value of the exponent 7 on EOS, it is 
found that the shape of the condensates at asymptotical 
times changes more dramatically. The non-canonical 
behaviour interpolating between 7=1/4 and 7 = for 
narrow resonances should be the easiest to identify. 



VI. CONCLUSIONS 



We have seen that an explicitly tunable two-channel 
model for cold Fermi gases with a narrow Feshbach reso- 
nance has a two-component order parameter. With each 
complex component having its own phase it is not sur- 
prising that, for the purposes of its EOS, the system can 
be described in a very transparent way by two coupled 
fluids. 



The system can be driven from one of Cooper pairs 
to one of tightly bound diatomic molecules by applying 
an external magnetic field. The single fluid limits of the 
deep BEC and BCS regimes, in which either the densities 
of Cooper pairs or molecules arc zero, are familiar. Our 
emphasis here has been on the very different single-fluid 
limit of the unitary regime, at which the scattering length 
diverges, in which neither density vanishes. For this we 
have shown that the EOS, with an exponent interpolating 
between 7 = and 7 = 1/4, will differ strongly from the 
canonical behaviour 7 = 2/3 of a conformal field theory, 
according as one, or both, of the effective range of the 
inter-atomic force and the length scale of the molecular 
field become large. This can be confirmed by measure- 
ments of the aspect ratio for elongated condensates. 



We conclude with an observation on formalism. Al- 
though the two-fluid description is very natural, the 
hydrodynamic equations can be reformulated as Gross- 
Pitaevskii (GP) equations for coupled complex GP fields 
*f = i/pFe** /y/2 and fij = = ^/p^e 18 * /y/2, one 
field for each density and phase. However, since phases 
arc field logarithms, this coupling is unmanageably loga- 
rithmic in the fields, in general. It is only when a single- 
fluid description is possible that the GP formalism is use- 
ful e.g. in the BEC regime, when *S? F oc A 0. 



A further case for which the GP formalism is useful, 
that we shall consider elsewhere [22|, is for the idealised 
situation in which [7 = 0, discussed in detail in [IcI ]. In 
this case we have a single-Quid description through the 
whole BEC-BCS regime or, equivalently, a single GP 
equation. The EOS has 7=1/4. 
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